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Abstract. Using the model of interaction between two-level systems (TLS) and strain 
fields introduced in Phys. Rev. B 75, 064202 (2007), we calculate the TLS-TLS interaction. 
We particularize our calculations to amorphous materials and analyze the dependence of 
the interaction Hamiltonian on the orientations of the two TLSs and of the direction that 
connects them. Assuming isotropic distribution of the TLSs orientations, we then averaged 
the interaction over one of the TLSs orientations. The resulting interaction depends in a 
simple way on the angle between the orientation of the other TLS and the line that connects 
the two TLSs. 
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1 INTRODUCTION 

Low temperature thermal properties of dielectric crystals differ markedly from 
those of dielectric amorphous solids. In crystals at low temperatures, the 
thermal properties are determined mainly by the acoustic phonons and their 
interaction with the defects; the system is well described by the Debye model 
and this renders the well known dependence of the specific heat (by T 
we denote the temperature) in a temperature range below 10 K. The heat 
conductivity of crystals is strongly dependent on the chemical composition of 
the crystal and on the nature and density of lattice defects. 

In contrast to the crystals, the specific heat of amorphous, glassy ma- 
terials below 1 K is proportional to the temperature, whereas the heat con- 
ductivity varies approximately like . Moreover-and this is maybe the most 
interesting characteristic of glassy materials-, some of their thermal properties 
(e.g. the heat conductivity) are (quasi) uniwersa/, i.e. they have a very weak 
dependence on the chemical composition or disorder of the solid. 

The thermal properties of glassy materials are described mainly in what is 
now called the standard tunneling model (STM) [HEj and have been attributed 
to the existence of dynamical defects in the solid. These dynamical defects 
are atoms or groups of atoms which tunnel between the two potential minima 
of quasi-symmetric two-well potentials. Because of the quasi-symmetry of the 
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potential landscape, the two lowest energy levels of the system may be very 
close to each-other and therefore may have a signicant contribution to the 
low temperature properties of the solid. If the potential barrier between the 
two wells is high enough, then the energy separation between the second and 
the third energy levels is much bigger than the energy separation between 
the lowest two energy levels and the system is well described by a two-level 
system (TLS), in a two-dimensional Hilbert space. In this Hilbert space, the 
Hamiltonian of a single TLS may be written as 



where az and ax are Pauli matrices. The parameters A and A in Eq. ([T]) are 
called the assymetry of the potential and the tunneling splitting, respectively. 

An elastic wave or a deformation of the solid perturbes the "free" TLS 
hamiltonian ([T]) by 



where in general we can write 5 = 2jijSij, where Sij is the strain tensor 
of the deformation field and "jij is a symmetric tensor that characterizes the 
TLS and its interaction with the deformation field-throughout this paper we 
assume summation over the repeted indices. While the tensor [7] is assumed 
to be the same for all TLSs, the parameters A and A vary form one TLS to 
the other, with a probability distribution 



where V is the volume of the solid and Pq is a constant. 

Although the STM describes well most of the properties of amorphous 
materials, it fails to explain certain "details", or some of its features need 
deeper understanding. One of the details that STM cannot explain in its 
simplest form is the fact that the heat conductivity is not exactly propor- 
tional to T^, but it is rather proportional to T^'®. Nevertheless, the most 
challanging problem related to amorphous materials is the explanation of the 
"universality" of their thermal properties. Why such a wide variety of mate- 
rials (polimers, oxide glasses, polycrystalline metals, etc.) have such similar 
properties and why the probability distribution P{A, A) has the simple form 
given by Eq. ([3]) for all of them? These questions are not yet answered, but 
there are hints that these effects are due to the interaction between the TLSs 
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P(A,A) = FPo/A, 



(3) 
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1.1 TLS-TLS INTERACTION IN THE STM 



The total hamiltonian of an amorphous sohd may be spht into three parts: 
the free phonons hamiltonian, the defects (i.e. the TLSs) hamiltonian, and 
the interaction hamiltonian between phonons and TLSs: 

H = Hph + Hdef + Hint, (4) 

where 

k,(T \ / 
Hdef = Y^^TLsim) (5b) 

m 

Hint = -YbijSij{m)az] (5c) 

m 

In Eq. ()5ap . the summation is taken over all the phonon modes, with k 
denoting the phonon wavenumber and a, the phonon's polarization; since we 
work with isotropic amorphous solids, the phonons will be considered as simple 
longitudinally or transversally polarised elastic waves (cr = / or t). Also in Eq. 
(|5a|) . M is the mass of the elementary "cell" of the medium, whereas pk.o- and 
Uk,(T are the momentum and displacement operators, respectively. 

In equations ()5bp and ([5c|) . the summations are taken over the TLSs in 
the system, therefore HTLsinT-) is the free Hamiltonian ([T|) of the m*^ TLS 
and Sij{m) is the strain field produced by the phonon system at the m^^ TLS. 

The interaction of the TLSs with the strain fields changes the ground 
state of the system, in the sense that the lowest energy state is not the one 
of zero strain, but one of non-zero strain. This non-zero strain that appears 
in the solid body at equilibrium, with the TLSs frozen in a configuration 
s(m), m = 1,2,..., induces long-range interaction between the TLSs. This 
interaction produces chaotic shifts of the energies of the TLSs, giving rise 
to spectral diffusion, which is the random change of the energy of a specific 
TLS due to the relaxation of the environment and is clearly manifested in low 
temperature hole-burning and phonon-echo experiments [3]. 

2 ANISOTROPIC TLS-TLS INTERACTION 

Minimizing the total energy dH) of the system one can obtain the interaction 
energy between the TLSs, 

Vl = Ui2aiz(T2z , (6) 
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where aiz and a2z are the Pauh matrices associated to the interacting TLSs, 
1 and 2, respectively. The couphng energy, U12, in a three-dimensional bulk 
material, has the expression (see for example [1]) 



U12 




(7) 




+6jknini + djiUiUk + SkiniUj) - IbniUjUkni] 



where by n we denote the direction from the TLS 1 to the TLS 2 and R is the 
distance between them; c; and q are the longitudinal and transversal sound 
velocities, respectively, p is the mass density of the solid, and 71,2 are the 
components of the tensors of coupling constants between the TLS 1 or 2 and a 
deformation strain field in the body. The second rank tensors will be denoted 
in general as [71,2]- From ([7]) follows that U12 is proportional to ii~^-like in 
the dipol-dipol interaction-and can be written as U12 = ui2/R^- Moreover, 
one can easily show that the average of U12 over the direction n, denoted as 
({7i2)n is zero for any 7i,2ij- Based on this, and since the expresion of U12 
is very complicated for general tensor elements, 7i^2ij~which are not known 
in detail anyway-, the usual assumption one does is that U12 is randomly 
distributed (i.e. there is no specific orientation dependence in the interaction) 
and a method similar to the random 3D Ising model is employed in analising 
the system [1]). 

In Ref. [5] we introduced a model in which a direction, call it t, is 
associated to each TLS and the TLS-strain field interaction depends on the 
relative orientation of t with respect to the strain field. Prom the components 
of t we construct the symmetric second rank tensor [T], with Tjj = titj and we 
write the [7] tensor in the general form, "fij = TkiRijki- Then the forth rank 
tensor [[i?]] is actually the tensor of coupling constants between the TLS and 
the deformation field and its general structure is determined by the symmetry 
properties of the solid in which the TLS is embedded. For an isotropic medium, 
[R\ is has only two independent parameters, C, and ^, similar to the Lame 
constants that enter in the construction of the elastic stiffness constants in 
elasticity: Rijki = [C ^ij^ki + ^' [^ik^ji + ^u^jk)]- Usually, for a direct connection 
to the STM, we take out a factor 7 = C' + 2^' and write Rijki = llC^ij^ki + 
ii^ik^ji+^ii^jk)], where C = C'/t and ^ = ^'/T! this implies = 1 Using the 

expressions for R and T, the [71,2] tensors become 7i,2ij = 7(C<^jj+2^ti, 24^1,2 j)- 
If we plug this into ([7]) and perform all the summations, we get the much 
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simpler expressions, 



U 



12 



^^{6e[(3C + 20cr' - 2(C + OcT'Eti ■ nf + (t2 • n)2] 

+ 12e2[5c,-2 - 4c^-2](ti . t2)(ti • n)(t2 • n) 

-60e2(q-2 - cr')(ti • n)2(t2 • n)^ - 4^2(3^-2 _ 2cr')(ti • t 

-4^(e + 2C)(q-2-c-2)-4Ceq'} 



(8) 



Using the relation + 2^ = 1, we elliminate C from all the expression above 
and get 



U 



12 



■nB?p 



2(l-Ocr1[(ti-fi)^ + (t2-ri 



+6e[5q 2 - Aci\ii ■ t2)(ti • n)(t2 • n) 



-2(2 -30 (q 



^)(ti-A) 



(t2 • n)' - 2C(3c, 
2(l-2e)q2} 



-2 



2cr')(ti-t2 



(9) 



In Eqs. dS]) and ©, L'^12 depends on the angles between ti, t2, and fi. The 
dependence is very complicated, but if we assume that in any volume element 
the TLSs are isotropically oriented, we can calculate the average interaction 
of TLS 1 with the TLSs located in a small volume 6V , at position ni? from 
TLS 1, by averaging over t2. This is simply done and we obtain 



12/ 



t2 



7'e(3 - 40 / 3__1 



(ti • fi) 



A_ 



(ti ■ n) 



(10) 



Using the relation 7^ = 4^7 /15, we can write the constant A as 



A = 




(11) 



We can observe from Eq. ()10p that the average of {Ui2)i^ over n gives zero, as 
one would expect. Moreover, the TLSs placed along the direction ti have-on 
averagean effect opposite to that of the TLSs placed in a plane perpendicular 
to ti-the plane that goes through the TLS 1. In other words, if {Ui2)i^ > 
at ti • fi = ±1, then {Ui2)^ < at ti • fi = 0, and viceversa. 

For concrete calculations of the interaction energy, we have to evaluate 
the parameter ^. This is a material dependent parameter and is determined 
by the ratio -fi/'jt [5], 

47? 15 



It 



40 

^-y + 32. 



(12) 
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Therefore, there is an ambiguity in determining ^. For each material, i.e. for 
every physically possible value of 7//7t, there are two solutions of Eq. (fT2]) . 
^1 and and up to now we have not been able to determine which one is 
physically relevant and which one is not; so in the calculations we have to 
consider both of them. But now from Eq. ()12p we notice that 



so for both values of ^, the modulus of A is the same and the ambiguity remains 
only in the sign: 



One of the most studied glassy materials is the amorphous silica (a-Si02 ). 
For this material 7; = 1.04 eV and 74 = 0.65 eV [6j, which give .^1 « 1.31 and 
^2 ~ 0.52. Using also the other physical parameters, q = 5.8 km/s, q = 3.8 
km/s, and p = 2200 kg/m^ [6], we obtain |^| 9.35 eVA^. 

As another example we take Epoxy. For this material, 7; = 0.35 eV, 
7t = 0.65 eV, Q = 3.25 km/s, p = 1200 kg/m'^ , while ct we approximate 
by the quasi-universal relation for glasses, q « q/1.65 [B]. Plugging all these 
parameters into the expression for A, we get |^| ~ 7.65 eVA^. 

3 Conclusions 

We used the model of Ref. [5] to calculate the TLS-TLS interaction in amor- 
phous glassy materials. In this model there are spatial directions associated 
to all the TLS. Therefore, the TLS-TLS interaction Hamiltonian that we ob- 
tained here depends on the orientations of the two TLBs that interact and also 
on the direction of the line that connects them. The angular dependence of 
the interaction allows for a more detailed treatment of the TLS-TLS interac- 
tion in glassy materials and a deeper understanding of the low temperature 
properties of these materials. 

Due to the isotropy of amorphous solids, one can assume that also the 
TLS orientations are isotropically distributed. We used this assumption to 
average over the orientations of one of the TLSs (say TLS 2) and we obtained a 
much simpler (average) interaction Hamiltonian, but which was still dependent 
on the angle between the orientation of the remaining TLS (the TLSl) and the 
line connecting TLSl to TLS2. Based on the knowledge of 7; and 7f, one can 
determine the interaction strength, up to its sign. Moreover, it was interesting 
to note that the interaction changes sign when going from the lateral direction 
(ti • n = 0) to forward or backward direction (ti • n = ±1). 
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